Lattice polaron formation: Effects of non-screened electron-phonon interaction 
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We explore the quasiparticle properties of lattice polarons on the basis of a quite general electron- 
phonon Hamiltonian with a long-range displacement-type of interaction. To treat the dynamical 
quantum phonons without significant loss of accuracy we adapt an exact Lanczos diagonalization 
method and compute various static and dynamical quantities, such as the electron-lattice correlation 
function, the polaron band dispersion, the effective polaron mass, the kinetic energy, the single- 
particle spectral function, and the optical conductivity, on finite one-dimensional lattices for a wide 
range of model parameters. We compare the results with those obtained for the standard Holstein 
model with short-range electron-phonon interaction only. 

PACS numbers: 71.38,+i, 63.20.Kr, 72.10.-d 
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I. INTRODUCTION 

The classical polaron problem [1] has received renewed 
attention on account of the observation of polaronic ef- 
fects in several important classes of materials, including 
high-temperature cuprate superconductors and colossal 
magneto-resistance manganites [2,3]. Remarkably even 
the much simpler case of free electrons interacting with 
optical phonons in ionic crystals is still not completely 
understood. From a theoretical point of view the chal- 
lenge is to describe the crossover from an only weakly 
dressed charge carrier to the strongly mass-enhanced, 
i.e., less mobile, polaronic quasiparticle with increasing 
electron-lattice coupling strength. Depending on the rel- 
ative importance of the short- or long-range electron- 
phonon (EP) coupling simplified models of Holstein [4] 
or Frohlich [5] type, respectively, have been studied over 
the last five decades. However, despite extensive ana- 
lytical work, in the physically most interesting crossover 
regime, up to now, the only reliable results came from nu- 
merical studies, such as finite-cluster exact diagonaliza- 
tions (ED) [6-11], (Quantum) Monte Carlo (QMC) sim- 
ulations [12,13], density-matrix renormalization-group 
(DMRG) approaches [14,15], and global-local [16] or vari- 
ational methods [17]. 

Recently the formation of small polarons was inves- 
tigated by Alexandrov and Kornilovitch [18] applying a 
new path-integral Monte-Carlo algorithm [19,20]. These 
authors introduced the following EP Hamiltonian 



Here c| and 6j TJ denote fermionic and bosonic annihi- 
lation [creation] operators, respectively. Restricting our- 
selves to the one-dimensional (ID) case, H describes an 
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(1) 



electron in a Wannier state on site j of an infinite chain 
which interacts with the vibrations of all ions of another 
chain via a "density-displacement" type long-range EP 
coupling 



(|i-j|3 + l)3/a 



(2) 



(cf. Fig. 1 of Ref. [18]). The distance |Z — j\ is measured 
in units of the lattice constant. In (1), xq = y / l/2Muio, 
kxq = y/EpUo , and the optical phonons, being polarized 
in the direction perpendicular to the chain, are approxi- 
mated as independent Einstein oscillators with bare fre- 
quency coq (h = 1). Physically, this model was proposed 
to mimic the interaction of doped holes with apical oxy- 
gens in the high-Tc's, e.g. in YBa 2 Cu306+ x , where one 
can assume that the coupling is not screened because of a 
low c-axis conductivity and high phonon frequency [18]. 
Methodically, model (1) represents an extension of the 
Frohlich model to a discrete ionic lattice or of the Hol- 
stein model including longer ranged EP interactions. In- 
deed, defining the polaron binding energy as 



1.27 e. 



(3) 



the Holstein model (HM) results by setting 
fiti) = K $j,i > 



(4) 



Therefore, the model (1) will be subsequently termed ex- 
tended Holstein model (EHM). In order to parametrize 
the EP coupling strength for both the HM and EHM we 
introduce two dimcnsionless EP coupling constants 

X = e p /2t, 

g 2 = s p /ujo (5) 
(in what follows we measure all energies in units of t). 
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So far, analytical and numerical investigations of the 
EHM have been mainly confined to the determination of 
the effective mass of the polaron, where it was found that 
the EHM polaron is much lighter than the small Holstcin 
polaron [18]. First results for the polaron band disper- 
sion and density of states were obtained quite recently, 
however, the QMC method of calculating the ground- 
state dispersion used by Kornilovitch [20] is limited to 
the case where the bandwidth is much smaller than the 
phonon frequency. 

In this paper we present an detailed comparative study 
of the Holstcin and extended Holstein models in order to 
discuss the effects of long-range EP forces on the lattice 
polaron formation. Using exact Lanczos diagonalization 
supplemented by a well-controlled phonon Hilbert space 
truncation method, we calculate for the first time spec- 
tral (optical) properties of the EHM polaron. As stated 
above such an numerical investigation is especially valu- 
able in the non-adiabatic intermediate-to-strong cou- 
pling transition region, where the electronic and phononic 
energy scales are not well separated, i.e., A ~ u> /t ~ 1. 
In the weak- and strong-coupling regimes the numerical 
work is supplemented by analytical approaches outlined 
in the Appendix. 



II. QUASIPARTICLE PROPERTIES OF LATTICE 
POLARONS 

A. Numerical Methods 

Before we discuss the various physical quantities 
let us briefly sketch our computational scheme. Di- 
agonalizing the coupled EP system (1) on hnite ID 
lattices with periodic boundary conditions (PBC), a 
general if -symmetrized state is given as \& K ) = 



E^oY:ILTcV\K;m,s), where S(m) = (N 1 + 
m)\/(N — l)!m!. K denotes the total momentum 
of the coupled EP system. Because the phonon 
Hilbert space has infinite dimension we apply a trunca- 
tion procedure restricting ourselves to phononic states 



\m,s) ph = n 



N-l 
1=0 



ph 



with at most M 



phonons, whereby m — X^o* n i — M, and nf e [0, m] 
(cf. Refs. [21]). The ground state 1$^ K=0 ) and all excited 
states \\P K ) contain components that correspond to m— 
phonon states in the tensorial product Hilbert space of 
electronic and phononic states. Accordingly, 



ml 2 



(M) 



S(m) 



=0 1 



(6) 



can be taken as a measure of the weight of the ra-phonon 
state in the K = ground state. In our ED analysis 
convergence is assumed to be achieved if the ground- 
state energy E n (M) is determined with a relative error 



less than 1CT 7 and \cff\ 2 {M) < lO -7 . Afterwards static 
correlation functions can be obtained easily by calculat- 
ing ground-state expectation values ($q(M)| . . . \)P (M)}. 
The numerical computation of dynamical properties, i.e. 
of spectral functions 



A°(u) = - lim -^m 

£^0+ IT 
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D-l 



u - H + E + is 
= ^ \($ n \Oi\# )\ 2 8[w-{E n -Eo)], (7) 

n=0 

is much more involved. Here O denotes the matrix repre- 
sentation of a certain operator O, and H is the very large 
sparse Hamilton matrix, acting in a Hilbert space with 
fixed momentum, which, for our problem, has a typical 
total dimension (D) of about 10 8 — 10 9 . Since it is im- 
possible to determine all the eigenvalues (E n ) and eigen- 
states (\&n)) of such a huge Hamilton matrix we com- 
bine kernel (Chebyshcv) polynomial expansion and max- 
imum entropy optimization in order to calculate A°(u>) 
in a well-controlled approximation (for more details see 
Refs. [21,22]). 



B. Electron lattice correlations 

In a first step we discuss the different nature of the 
polaronic states in the HM and EHM in terms of static 
correlation functions (rnqi) between the electron position 
[i = 0] and the oscillator displacement [qi cx (b\ + b t )] at 
site I, 



Xo.i = (clc {bl +l + b 0+l ))/N . 



(8) 



X i indicates the strength of the electron induced lattice 
distortion at i = and its spatial extent [14,23,17], where 
Af = X)/( c o c o(^o+/ + ^o+;)) ^ s a normalization constant 
(note that Af = 2(e p /ujo)(c' c ) holds for the HM). 

Figure 1 shows the (static) electron-lattice correla- 
tion function (8) in the weak- (a) and intermediate-to- 
strong- (b) EP coupling regimes, where we have cho- 
sen an intermediate phonon frequency (loq — 1) in order 
to include non-adiabatic effects. Clearly for the quan- 
tum phonon model (1) the EP interaction gives rise to a 
"dressing" of the charge carrier at any finite A, g 2 . If the 
EP coupling is weak, however, the amplitude of x l is 
small V I (in particular smaller than the quantum-lattice 
(zero-point) fluctuations), that means the lattice defor- 
mation could not trap the charge carrier and a so-called 
"large" polaron (LP) is formed in both the Holstein and 
extended Holstcin models. Obviously, the situation is en- 
tirely different in the strong-coupling region. For the HM 
the EP correlations are almost local indicating the for- 
mation of a "small" polaron (SP). On the other hand, as 
a result of the non-screened EP interaction, in the EHM 
the deformation is spread over many lattice sites, i.e., we 
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found again a LP. It is worthwhile to point out, however, 
that the electron and the phonon cloud are tightly bound. 
That means the LP of the EHM as a whole behaves as a 
well-defined polaronic quasiparticlc (cf. Sec. II C) and, 
in our opinion, it is not possible to discuss the size of 
the electronic wave function and the size of the lattice 
distortion separately [18]. 

In the insets of Fig. 1 we show the differences be- 
tween the phonon distribution functions in the weak- 
and strong-coupling cases, where the ground state is ba- 
sically a zero-phonon- and multi-phonon state, respec- 
tively. With regard to the discussion of the effective 
mass in Sec. II E we would like to annotate here, that 
at small (large) A the EHM polaron contains more (less) 
phonons in its phonon cloud than the HM polaron. Of 
course, in the extreme strong-coupling limit, the usual 
Poisson distribution with parameter g 2 results, demon- 
strating that adjusting the parameters of both models 
according to (3)-(5) is correct. 

C. Single— particle spectral function 



which is only slightly reduced from the free electron 
value, the wave-function renormalization factor Z^ =1T is 
almost zero. The WCT shows that the state with K = ir, 
being energetically separated by wo from the ground- 
state energy, is predominantly a phononic state. At 
strong EP coupling the polaronic band is characterized 
by Zk <C 1 VK, indicating a strong mixing of electronic 
and phononic degrees of freedom. Calculating the pola- 
ronic quasiparticlc weight factor within the framework of 
the strong-coupling theory (SCT) developed in Appendix 
B (Eq. (44)) gives Z K=0 -v&\ues which are by a factor of 
3 too small as compared to the exact data of Fig. 2 (b) 
and (d). The differences mainly arise because these pa- 
rameters correspond rather to the intermediate-to-strong 
than to the extreme strong-coupling regime. The qual- 
itative behavior of the single-particle spectral function 
Ak{E), however, is correctly reproduced by SCT (43), 
which yields, above the quasiparticle pole, a sequence 
of excitations separated by ujq. Apparently the spectral 
weight of this incoherent part increases with increasing 
EP interaction strength. 



Next, in order to examine dynamical quasiparticle- 
properties of the HM and EHM polarons, we have eval- 
uated the wave-vector resolved spectral density function 

A K {E) = ]T | ($ niK | 4 | 0) | 2 6(E - E n , K ) . (9) 

n 

The results are presented in Figure 2. To visualize the 
spectral weights of the various excitations, the integrated 
density of states, 

N K {E) = / cLE'Ak(E'), (10) 

J — oo 

is also displayed. The weight of the first delta-function 
peak in each If -sector gives the wave-function renormal- 
ization factor [24] 

Z k = \(%,k\J k \0)\ 2 , (11) 

where \&o,k) denotes the single-polaron state with mo- 
mentum K being lowest in energy. Z K=0 is usually 
termed "quasiparticlc-weight factor" . Since the total in- 
tegrated area under the entire spectra is unity, the renor- 
malization factor is less than unity and, in particular, 
Zk=o is a measure how much the polaronic quasiparticle 
"deviates" from the free electron (Zk=o = 1)- In ac- 
cordance with the discussion in the preceding section, in 
the weakly-interacting EHM we found a stronger dress- 
ing of the electron by phonons than in the HM, i.e., a 
larger renormalization Z^^m < Z^^o- This finding is 
corroborated by the the weak-coupling theory (WCT) 
outlined in Appendix A. Table I demonstrates the good 
agreement of the theoretical approach, working for the 
infinite system, and finite -cluster diagonalizations, pro- 
vided that both A and g 2 are small. Contrary to Zk=o, 



D. Polaron band structure 

Now the so-called "coherent" band dispersion, Ek, 
can be derived from the first peak of Ak (E). Figure 3 
compares the dispersion of the energy bands for the EHM 
at different EP interactions, corresponding to the weak-, 
intermediate- and strong-coupling case, where, by going 
from (a) to (c), uo is increased modelling the adiabatic, 
intermediate, and anti-adiabatic regimes. 

Starting with the adiabatic weak-coupling case 
(Fig. 3 a), we found that the band structure is nearly 
unaffected at small momentum, i.e. in the vicinity of 
the band center. In this region I^q^) is quasi a zero 
phonon state [10,26]. A different behavior is observed 
near the zone boundary. Here the band structure is flat- 
tened. Such a "flattening" has been found for the HM 
as well and can be attributed to the intersection of the 
dispersionsless optical phonon branch with the bare elec- 
tronic cosine band [25,9,10,26] (cf. also Fig. 4). The 
weak-coupling calculation of Ek (Appendix A; Eq. (22)) 
reflects this behavior. For the HM, e.g., the correction to 
£if is given by the integral X^(K, 0), which is non-zero 
only for (Ek — ^o) 2 > 4i 2 . The latter condition yields 
a threshold K* , at which the solution of (22) jumps to 
the bare band dispersion £k- For K > K* , the first ex- 
citation in A k (E) is related to a one-phonon absorption 
process and as a result the WCT approximation for Ek 
breaks down. Thus, above K* , the physical solution is 
given to lowest order by the dashed line at Eq + ujq. 

The flattening considerably weakens and ultimately 
vanishes if the EP coupling A increases. This tendency is 
especially pronounced for the EHM in the non-adiabatic 
regime. As can be seen from Figs. 3 (b) and (c), in the 
strong-coupling non-adiabatic regime our SCT yields ex- 
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ccllent results. Most notably we do not observe the same 
drastic polaronic band collapse as in the HM [10,11], i.e., 
in the EHM the coherent bandwidth AE = E n - E 
becomes much less renormalized by the EP interaction 
(e.g., for the HM with A = 5.0 and uoq = 3.0, we found 
AE = 0.15319). 



E. Effective mass 

Another important question is the change in the po- 
laron effective mass induced by the EP coupling. In 
general it is difficult to compute the mass enhancement, 
which is defined as an inverse second derivative of the 
band energy with respect to quasi- momentum at the 
band minimum, m*/m oc [d 2 E(K)/dK 2 \K=o]~ 1 , using 
finite-lattice diagonalizations, because E(K) is known at 
multiples of 2ir/N only rather than at any K, making the 
limiting procedure if —> ill posed. On the other hand, 
the mass enhancement factor m* /m is also related to the 
quasiparticle weight factor Zk=q [24]. For the HM, we 
are able to prove the relation 



m/m 



HM 



7 HM 



(12) 



in the weak-coupling limit (see Appendix A, Eqs. (23)- 
(25)), i.e., at A <C 1 the polaron effective mass can be 
read off from the first step in the integrated spectral 
weight function depicted in Fig. 2. Plotting Z^^ as 
a function of A in Fig. 5 and comparing the effective 
mass determined in this way with the QMC masses, ob- 
tained from m/m* — d 2 Ex / OK 2 without any systematic 
finite-size errors [27,18,20], the perhaps surprise finding 
is that Eq. (12) holds for the whole coupling region. That 
means, in the Holstcin model, we can determine the effec- 
tive polaron mass simply by calculating the quasiparticle 
weight factor. In previous ED studies of the Holstcin 
polaron problem this fact has been ignored so far. 

Unfortunately no such simple relation exists for the 
EHM. This is shown more explicitly in the Appendix, 
where approximative expressions for m/m* were derived 
in the weak- and strong-coupling limits. Note that our 
analytical weak- and strong-coupling approaches con- 
firm the unexpected non-monotonic dependence on A of 
m *EHM l m *HM i which was found numerically by Alexan- 
drov and Kornilovitch [18] (see inset). In the light of 
the results presented in the previous sections, it becomes 
clear that at small EP couplings the EHM LP has to 
drag a larger phonon cloud coherently through the lat- 
tice than the HM LP and therefore acquires a larger effec- 
tive mass. Further numerical data show that this effect 
becomes negligible in the weak-coupling anti-adiabatic 
regime, where the phonons can follow the electron in- 
stantaneously (cf. also Fig. 3 of Ref. [18]). As a matter 
of course, in the strong-coupling limit the EHM LP is 
much lighter than the HM SP due to the weaker band 
rcnormalization caused by the extended form of the lat- 
tice distortion. 



F. Optical conductivity 

In this section we compare the optical response of HM 
and EHM polarons. The real part of the optical conduc- 
tivity, 



Re<7(w) = V5{lo) +<r reg (uj) , 



(13) 



can be decomposed into the Drude term (oc T>) at lo = 
and a regular contribution for lo > 0, which, in linear re- 
sponse theory, for the (extended) Holstein model is given 

by 



reg 



ao \^ \(%\itEj(c]c j+1 - ct +lCj )|^J| s 
N 



E 

n#0 



E n — Eq 

x 5[u> - (E n - E Q ) 



(14) 



with do = 7re 2 (T = 0, K = 0-sector). Again we intro- 
duce a w-integrated spectral weight function, 



put 

S reg {oj) = / duj'<j reB {J) , 

JO 



(15) 



in order to visualize the intensity of the various excita- 
tions more clearly. 

Figure 6 shows the optical conductivity obtained at 
loq = 1.0 for the HM and EHM on a eight-site lattice 
using PBC. For the ID HM the optical absorption spec- 
trum has been discussed in detail in previous work [23]. 
If the energy to excite one phonon lies inside the bare 
tight- binding band we found, at weak EP coupling, the 
first transitions by adding phonons with opposite mo- 
mentum to only slightly renormalized electronic states 
(in order to reach the K = sector of the ground state). 
However, since the ground state is approximately a zero- 
phonon state (cf. inset of Fig. 1 (a)), the spectral weight 
of optical transitions involving larger number of phonons 
is reduced drastically. Of course, the absorption thresh- 
old is loq for the infinite system; the shift observed in 
Fig. 6 (a) simply results from the discrete if-mesh of 
our finite system. The situation changes by increas- 
ing the EP coupling when in the HM the SP formation 
takes place. Now the phonon distribution function in the 
ground state is broadened and in the optical response the 
overlap with excited multi-phonon states is enlarged. As 
a result the famous SP absorption maximum develops 
around lo ~ 4A = 2e p for large couplings. 

Let us now discuss the optical response in the EHM. 
Of course, there is little change in the weak-coupling 
region. At large EP coupling, however, the optical ab- 
sorption points toward a completely different nature of 
the polaronic states in the Holstein and extended Hol- 
stein models. The EHM polaron clearly shows all the LP 
signatures but compared to the weak-coupling case and, 
what is more important, also compared to the HM SP, 
the optical conductivity is strongly enhanced due to by 
multi-phonon absorptions processes. The physical rea- 
son lies in the non-screened EP interaction leading to the 
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form of the lattice distortion depicted in Fig 1 (b) . Tak- 
ing into account the internal structure of the EHM LP it 
is obvious that the lattice distortion undergoes less rel- 
ative changes when the charge carrier hops incoherently 
to neighboring sites accompanied by phonon absorption 
or emission [18]. 

G. Kinetic energy 

Integrating (13) with respect to u>, the familiar f-sum 
rule 

E km S tot V S re ° 

/ (To ^CTo (To 

can be derived, where S reg = S re9 (oo). Eq. (16) relates 
optical response and kinetic energy. E kin measures the 
mobility of the charge carrier. 

To elucidate the different nature of HM and EHM po- 
larons in more detail, in Fig. 7 we have displayed the ki- 
netic energy (oc S tot ), renormalized to its value at A = 0, 
together with S reg . Since the kinetic energy contains 
contributions from both "coherent" (oc T>) and "incoher- 
ent" (oc S rea ) hopping processes, the Drude part can be 
directly read off from the difference between the filled 
and open symbols at fixed A. In agreement with pre- 
vious numerical results, in the HM we found a contin- 
uous transition from a LP to a less mobile SP as the 
EP interacting increases [12,23,14]. The decrease of S tot 
in the crossover region, being much more pronounced in 
the adiabatic regime [23] (as well as for higher dimen- 
sions [12,11]), is driven by the sharp drop of the Drude 
weight [28]. By contrast, in the EHM the kinetic energy 
decreases very gradually with increasing A and we ob- 
serve a substantial Drude contribution even at large EP 
couplings. This is in accord with the moderate renor- 
malization of the polaronic bandwidth and the minor ef- 
fective mass (cf. Sec. II D and Sec. E). In addition, as 
already stressed in Sec. II E, the optical absorption due 
to inelastic scattering processes, described by the regular 
part of the optical conductivity, gives a large contribu- 
tion. This can be easily understood within second-order 
perturbation theory (note that one has to go beyond the 
lowest order of approximation to obtain reliable results 
for the kinetic energy [29,9,11]; cf. Appendix B): during 
a second-order hopping process of the EHM LP the lat- 
tice distortion of the intermediate state with the charge 
carrier on a nearest-neighbor site of the initial site fits 
much better to the polaronic quasiparticlc than in the 
case of the HM SP. The difference between the numerical 
and theoretical results at larger EP couplings originates 
from the neglect of longer-ranged hopping processes in 
our theoretical approach. Of course, such transitions are 
much more important in the EHM. In addition, let us 
emphasize that S^lj M > holds also in the weak- 

coupling limit (see inset). That means the stronger re- 
duction of the coherent Drude part, corresponding to the 



stronger mass enhancement in the weakly-coupled EHM 
( m *E h m / m *H m > 1' c f- inset Fig. 5), is overcompensated 
by the incoherent part. 

III. CONCLUSIONS 

In summary, we have performed an extensive compara- 
tive numerical study of polaron formation in the Holstcin 
and extended Holstein models, supplemented by a theo- 
retical analysis of the weak- and strong coupling limits. 
The emphasis was on the new effects induced by a non- 
screened electron-phonon interaction. The main charac- 
teristics of the new polaronic state formed in the EHM 
are the following. 

(i) By its nature the EHM polaron is a large polaron 
in the whole EP coupling region. That is the lat- 
tice distortion is spread over large distances even if 
the EP is extremely strong, In this regime a small 
polaron is formed in the Holstein model. 

(ii) For strong EP interactions the EHM polaron prop- 
agates in a relatively weakly renormalized band 
as compared to the HM. Accordingly the effective 
mass of the large EHM polaron is much smaller 
than that of the small Holstein polaron with the 
same polaron binding energy. 

(iii) A surprise finding is that the effective mass of the 
EHM polaron, describing a "coherent" band mo- 
tion, is larger than the effective mass of the HM 
polaron at weak EP couplings, in particular in the 
adiabatic case. We have seen that this effect can 
be attributed to the larger number of phonons the 
charge carrier has to drag through the lattice if the 
weak EP interaction is non-screened. 

(iv) From the calculation of the if-resolved single par- 
ticle spectral function a wave-vector renormaliza- 
tion factor Zk can be extracted, which indicates, 
in accordance with (i) and (ii), at weak (strong) 
EP couplings a stronger (weaker) renormalization 
of the band states in the EHM than in the HM. 

(v) While in the HM the inverse polaron effective mass 
is directly given by the quasiparticle weight fac- 
tor, Zk=o, the relation is more complicated for the 
EHM. In the weak-coupling limit this has been cor- 
roborated analytically. 

(vi) The EHM polaron band dispersion is non-cosine 
for model parameters corresponding to the (adia- 
batic) weak-to-intermediate coupling regime. In 
particular, in the weakly-coupled EHM, a flatten- 
ing of the band structure at the zone boundary is 
definitely observed just as in the HM, but the effect 
is much less pronounced. Furthermore the flatten- 
ing rapidly vanishes with increasing EP coupling 
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strength and phonon frequency. In the strong- 
coupling limit, the EHM exhibits a free-particle 
like band dispersion with a bandwidth which, al- 
though renormalizcd, is approximately one or two 
orders of magnitudes larger than in the HM. 

(vii) While in the HM the transition from large to small 
polarons is accompanied by significant changes in 
the optical response, the optical absorption in the 
EHM shows large polaron characteristics for all EP 
interaction strengths. Most notably the extended 
form of the the lattice distortion in the EHM gives 
rise to a large amount of "incoherent" hopping pro- 
cesses contributing to the regular part of the optical 
conductivity. As a result the regular contribution 
to the f-sum rule is always bigger than in the HM. 

(viii) If one takes the averaged kinetic energy as a mea- 
sure for the mobility of a charge carrier, the EHM 
polaron is more mobile than the HM polaron, in- 
dependently of the magnitude of the EP coupling 
strength. In particular the dramatic kinetic energy 
loss during the self-trapping transition of the Hol- 
stein small polaron is absent in the EHM. On the 
contrary, one observes a very gradual decrease of 
the kinetic energy with increasing EP interaction 
and a substantial Drude contribution even for large 
EP couplings. 

Finally we would like to stress that the above properties 
of the EHM large polaron are generic and not an arti- 
fact of our ID system. The relative mass enhancement 
m EHM/ m *HM a t weak EP interactions, e.g., seems to be 
even more pronounced in 2D [18]. 
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APPENDIX: ANALYTICAL APPROACHES 

Following the previous consideration of the HM [11], 
the Hamiltonian for treating the EHM at any A may be 
written as 

H = -v E c h E c o'A> c i + ^ E(^ + \) > 

3 3,3' I 

where r\ is a c-number and Cjij are generally functions of 
the phonon operators b\, b v Using the formalism of gen- 
eralized Matsubara Green's functions, the polaron self- 
energy S corresponding to (17) in the second step of 
iteration is given by 



E(jm; j 2 r 2 ) - -{C n32 )5{ri - r 2 ) + ]T Q(j'r 1 ;j"T 2 ) x 

j'j" 

x [(T T C, U , (tOC,^ (r 2 )) - {C hj ,){C r , j2 )] , (18) 

where Q(j' T\\j" Ti) means the polaron Green's function 
in the first approximation (see Refs. [30-32] for details). 
In the subsequent calculations, the latter equation will 
be taken as a starting point for the treatment of EHM in 
the weak (A <C 1)- and strong (A 3> l)-coupling regimes; 
in addition, the low-temperature approximation (LTA) 
defined by /3wo 3> 1, and the small carrier-concentration 
limit (x — ► 0) will be assumed. 



A. Weak— coupling regime 

In this case, the functions Cj>j are defined by the first 
and third terms of the Hamiltonian (1), and r\ is put 
equal to the chemical potential fi. Applying the Fourier 
transformation to both sides of (18) and carrying out 
the standard summation over the phonon Matsubara fre- 
quencies [24,11], the polaron self-energy 

oo 

E K (u>) = U+^o £p( d ) cos(Kd) (19) 

d— — oo 



■ + [1 - LU - £ K , 



is obtained after analytical continuation to the real fre- 
quencies uj under assumption of LTA and small x. In 
writing (19), the definition 



(20) 



£k = —2tcosK, and w = w + i0 + were used. The HM 
limit results from (19) by setting 



£ P {d) — £pS d Q . 



(21) 



In view of (19), the polaron band energies Ek are so- 
lutions of the following equation 

oc 

Ek=Zk+wo E £ P (d) cos(Kd)lW(K,d), (22) 

d— — oo 

and the renormalization factor of the spectral func- 
tion [24], Z K = [l--^ReE K (cj)]~t EK _^ is deter- 
mined by 



Z K X = 1 + uj a i p(d) cos(Kd)lM(K,d) , (23) 

d— — oo 

where 
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The relation between the effective polaron mass to* 
and the bare electron mass m (being equal to (2t)^ 1 in 
ID) is deduced according to to/to* = [8Ek /deK]\s K ^o 
with e K = tK 2 [24]: 



— = Z 

TO* 



K=0 



d= — oo 



(25) 



Note that only for the Holstein model (d — 0) the effec- 
tive mass is given by the inverse spectral weight factor 
TO HM/ m = (Zr=o) ■ Using (22), the polaron kinetic 
energy results from the relation E^in = tdtE^K^a as 



-2tZ 



1+UJ £ 



(26) 



n dK' cos if' cos(K'd) 
„ 2tt 



B. Strong— coupling regime 

In order to generalize the strong-coupling approach, 
developed in [11] for the HM, to the EHM case, the long- 
range interaction of Hamiltonian (1) is eliminated by a 
non-local Lang-Firsov transformation 



u = U c ^{^T,f^H c M -h) 



u> 



(27) 



Clearly, the theory based on (27) turns into the theory 
for the HM if the condition (4) is inserted. The canonical 
transformation (27) applied to (1) leads to the polaron 
binding energy (3) and to the emergence of the multi- 
phonon processes connected with the electron hopping 
from the site j to the nearest-neighbor sites j + h (h 
being the elementary translation in units of the lattice 
constant). 

Treating the dynamical EP interaction of the trans- 
formed Hamiltonian by means of the formalism outlined 
in [11], the polaron self-energy represented in the space 
of Brillouin-zone if -vectors and Matsubara frequencies 
iujv = i(2v + l)n/(3 is obtained to the second order as 



3', 3" 
32-31 



x l y e iK'w-i") \ f dT s»»-»^ (28) 

N k^c 13 -* 

x i 2 (feyW^ j2 (0)) - (1> jl , j >)(<l> j », j2 ))GK>(iu C ) ■ 



The multi-phonon operator 



$ jlj+h = cx P j ^ £ {-M + h)- Mj)) (b] k) 



(29) 



occurring in the transformed hopping term of (17) in- 
duces the polaron band narrowing in the first order, 
namely, 

i=t($ j>j+h ) = icxp{-g 2 A(l)coth(/?w /2)} , (30) 
where 



A(l) 



i -£/«(!)/ £/ ; 2 (o) 



(31) 



In the LTA f3uj > 1, we have t ~ t exp{-g 2 A(l)}, £ K 
—2tcosK, and 



-i") 



3' ,3" 
32 -31 



P£ 



X s 1 



2sujq 



^ s\ (3 (su ) 2 + (uj ( - u v ) 2 
The parameter k = x{d, d\,d 2 ) is given by 

x = 4 £ M°) \M d ) + M d + <h-di) 

^0 J 

-fl(d-d 1 )-f l (d+d 2 ) . 



G K '(icoc)- (32) 



(33) 



Here d = j 2 - ji, d\ = j' - j\ and d 2 = j" - j 2 are 
elementary translations. The summation over ( is now 
evaluated using the lowest-order approximation for the 
Green's function Qx(i^c) corresponding to a quasiparti- 
cle energy spectrum of the form £k — e p . After carrying 
out this summation in LTA and for small carrier concen- 
tration, the wave-vector- and frequency-dependent po- 
laron self-energy is obtained from (32) by the analytical 
continuation iuj v -^was 



d.di.do K' 



(34) 



£ 



-iK'(d+d 2 -d 1 ) _ 



SI 



uj - £k' + £ P + M _ SLU o 



If we neglect £k' on the right hand side of (34), in the 
strong-coupling limit, the polaron band dispersion Ek 
can be easily determined from 



Ek = h - 2t 2 £ 

+^2K£ 



(2 5 2 A(1)) S 



1 



s\ suj q - E K 

S>1 

[g 2 A(2)Y 1 



(35) 



S>1 



si suj a — Ek 



where 
A(2) 



, 2^/ ; (Q)/Ki)-E ; /KQ)/K2) 



(36) 
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On the basis of (35), the polaron mass enhancement 
can be calculated from 



m/m* 



-9 2 A(1) 



[dE K /de K ]\e K ^o ■ 



(37) 



where Ek = tK 2 . Substituting, for K — > 0, £k — —2t + 
Ek and i&K = — 2i(t — 2ek), we find 



m 
m* 

with 

Z Ek 



(38) 



1 + 4te g 2 (A(2)-A(l)) 



SWo - ^0 / S 2 A (2) 



i + 2r 



+e 



,^(A(2)-2A(1)) 



(SW - Sk) 2 / 2 g2A(l) 
1 



(39) 



(sw - E K ) 2 / „2 



cos2A 



9 2 A(2) 



Here (. . .) e denotes the average over s > 1 with respect 
to the Poisson distribution with the parameter g. 

Accordingly, the polaron kinetic energy takes the form 



Ekin — ~^Z E } 



e - 3 2 A(l) + 2t 



+2te g2(A(2) - 2A(1)) 



1 



suj -E / 2 S 2 A (i) 

sujq - E / g2A{2) 



(40) 



Finally, in order to discuss qualitatively the behavior of 
the single-particle spectral function reported in Sec. II C, 
we calculate Ak{u>) in the strong-coupling limit. Apply- 
ing the transformation (27), the electron operators are 
transformed into 



transformation of Eq. (42). As a result, the spectral func- 
tion 



A K (u) = -2ImC?f (w) 

= e- fl2 Z^2n5(u - [E K -e p ~n]) 



(43) 



EE^Ei(^) -<a--a-v] 



A ^ S! 

d s>l AT' 



w 



:Z E ^,27r5(w - [E K i + suj -£ P - lA) 



is determined by the spectral functions Ak>{u — ru>o) 
with r > corresponding to Gj^'i^ — rw o), i- c -, using 
the self-energy Ejfi(w') derived in this section. The en- 
ergies on the r.h.s. of (43) are solutions of (35) and 
the factors Ze k , are given by (39). The first term on 
the r.h.s. of (43) describes the quasiparticle of momen- 
tum K and the second term, being a sum over the entire 
polaron band, corresponds to the incoherent part of the 
spectral function. Inserting the condition (21) into (43), 
the spectral function of the Holstein model is obtained, 
which differs by the self-consistently determined Ek' and 
Ze k , from the result of Alexandrov and Ranningcr [34] . 

In particular, at K = 0, r = 0, the quasiparticle weight 
factor results as 



J E 



(44) 



Here it is necessary to point out the different level of 
approximation we used in deriving Eqs. (38) and (44). 
Therefore it is not possible to verify the HM relation (12) 
by the above strong-coupling calculation of Zx=o- The 
leading exponential dependence of (38) and (44), how- 
ever, is found to be the same (not the same) for the HM 
(EHM), in good agreement with the numerical results of 
Figs. 2 and 5. 



=(t) 



cxp { ( 



o^EM')^- 6 *) 



(41) 



Consequently the spectral function is determined by the 
imaginary part of the retarded Green's function G R (K, to) 
of the operators c K , c* K . Owing to the relation between 
the time-ordered products of operators and the retarded 
Green's functions of the same operators [33] we consider 
the time ordered product (T t c.K(t)c* K (Q)) and perform a 
decoupled average over the phonon and charge variables. 
In the limit T-^0we get 



<T t c*(t)4(0)> 



(T tCjf (t)4(0)) 



N 

EEE 

K' d s>l 



1_ / £p(d) \ —isu a t i(K — K 



'W(T t CK>(t)cU0))]- (42) 



The relation between the imaginary parts of retarded 
Green's functions G r {K,oj) and Q r (K 7 uj) of operators 
Ck and ck, respectively, is obtained from the Fourier 
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TABLE I. Quasiparticle weight Zk=o obtained from ED 
(N — 8, M — 24) and within WCT according to Eq. (23). 
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0.1, g 2 = 0.2 


A = 0.5, g 2 


= 1/3 
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WCT 


HM 
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0.946 


0.893 


0.848 


EHM 


0.918 


0.893 


0.857 


0.781 



FIGURE CAPTIONS 

Fig. 1: Electron-lattice correlations in the weak- (a) 
and strong-coupling (b) cases. ED results arc obtained 
for a finite chain with N = 8 sites and at most 24 
phonons. The insets show the weight of the m-phonon 
state in the ground state. 

Fig. 2: Single-particle spectral function Ak(E) (thin 
lines) and partial integrated density of states Nj^(E) 
(thick lines) for the ID HM (a-b) and EHM (c-d) with 
uj = 1.0 (N = 8, M = 24). Solid and dot-dashed lines 
belong to states with total momentum K = and K = n, 
respectively. 

Fig. 3: Band dispersion of the ID extended Holstein 
model for low (a), intermediate (b), and high (c) phonon 
frequencies. Exact data are extracted from finite-lattice 
diagonalizations with N = 8 and 10 sites. In the weak- 
and strong-coupling regimes ED results are compared 
with the theoretical predictions. 

Fig. 4: Flattening of the polaron band dispersion in 
the ID Holstein model (weak-coupling case). Long and 
short-dashed curves denotes the solution of Eq. (22) and 
the bare phonon frequency uj = 1.0, respectively. The 
agreement of ED and WCT gets better as A decreases. 

Fig. 5: Inverse effective polaron mass for the ID (ex- 
tended) Holstein model at u>o = 1-0. The QMC data are 
taken from Ref. [18]. The ratio of the effective masses of 
the EHM and HM polarons are displaced in the inset. 

Fig. 6: Optical absorption in the ID Holstein and 
extended Holstein models. The regular part of the con- 
ductivity a reg (thin lines) and integrated spectral weight 
greg ^ mc ]j h ne s) are shown in the weak- (a) and strong- 
coupling (b) regimes. 

Fig. 7: Rcnormalizcd kinetic energy (S tot ) and contri- 
bution of a re9 to the f-sum rule (S re9 ) as a function of 
EP coupling (A) at ui n = 1.0. 
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